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Abstract 

For twenty years research into the anomalies in the HF spectra was going in a 
wrong direction in fighting the Bohr—Weisskopf effect. As way out, we propose 
the model-independent way, which enables the nuclear radii and their moments to 
be obtained from the hyper-fine splitting. The way is based on analogy of HFS to 
internal conversion coefficients, and the Bohr—Weisskopf effect — to the anomalies 
in the internal conversion coefficients. It is shown that the parameters which can 
be extracted from the data are the even nuclear moments of the magnetization 
distribution. The radii i?2 and (for the first time) R4 are obtained in this way 
by analysis of the experimental HFS for the H- and Li-like ions of 209 Bi. The 

critical prediction is made concerning the HFS for the 2 pi state. The moments 

2 

may be determined in this way only if the higher QED effects are properly taken 
into account. Therefore, this set of the parameters form a basis of a strict QED 
test. Experimental recommendations are given, aimed at retrieving data on the 
HFS values for a set of a few-electron configurations of various atoms. 


1 



1 Introduction 


Hyperfine structure e.g. of the D\ line of the 133 Cs atom or that of the Is — 2s tran¬ 
sition of the hydrogen atom play an important role in construction of the atomic 
clock (e.g. Cl)- Experiments with the hyperfine structure allowed one to obtain 
more accurate value of the hyperfine structure constant a — one of the fundamen¬ 
tal constants. Worthy of notion is also the project of composing a reference point 
of frequency, founded on the few-electronvolt isomeric state of the 229 Th nuclide [2], 
with unprecedented systematic shift suppression, allowing for the atomic clock per¬ 
formance with a total fractional inaccuracy approaching 10~ 19 - 10 -21 . Hyperfine 
structure is of especially great importance for the proper design of such clock [3j. 
At the same time, proper account of the hyperfine structure may radically change 
the estimated lifetime of the nuclear isomeric states [4j. 

From the viewpoint of theory, the hyperfine structure arises due to interaction of 
the electron magnetic moments with the nuclear magnetic moments. The interac¬ 
tion leads to a hyperfine splitting (HFS). That also depends on the electron density 
near the nucleus and other properties of the electron shell. In the first approxima¬ 
tion, both nuclear and atomic factors act independently, which manifests itself in 
the factorization of the nuclear and electron parameters in the related formulae. 
However, precision measurements showed that the factorization is violated for the 
reason of the finite nuclear size. First of all, this is the Bohr-Weisskopf effect: finite 
distribution of the magnetic currents inside the nucleus. Study of the magnetic 
anomalies offers a way of pursuing the change of the effective nuclear radii within 
isotopic chains with addition of neutrons to the nucleus, which information is of 
primordial interest. That explains why a specific problem of a magnetic anomaly 
of the hyperfine splitting of the atomic structure became very topical for the past 
time B El 0 i- Nuclear-optical methods became very important in the research 
into the properties of rare or radioactive nuclides. 

From the theoretical point of view, description of the hyperfine splitting is a 
challenging problem, as it needs a calculation of the atomic wavefunctions with the 
accuracy of ~ 10 -4 [5j[9]. A convenient way of comparison of the data to theoretical 
calculations is offered by investigation of the HFS in few-electron atoms: H-, Li- 
or B-like ions. In these cases, the electronic wavefunctions can be calculated more 
reliably. Moreover, comparison to the theory can be used as a test for our capability 
of description of both the electronic structure of the atoms, as well as QED. In 
this aspect, the problem of the Bohr-Weisskopf effect in appearance of the magnetic 
anomaly is put forward in the first place. Existing method of its description through 
construction of the specific differences is however far from being appropriate and 
sufficient for this role. It is erroneous in principle, as it is clear from the internal 
conversion theory and is shown in the present paper. 
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An adequate way of solving this problem was proposed in refs. mmm - An 
analogy was noted between the internal conversion coefficients (ICC) and the HFS. 
It was shown that the HFS can be calculated as the limiting case of the ICC at 
the transition energy ui —>• 0. Test calculations were made for the H-like ions of 
209 Bi. Our present purpose is to show that the application of the classical theory 
of anomalous internal conversion H3 03] gives us the proper tool for adequate 
description of the Bohr-Weisskopf effect, also offering the way of extracting the 
unique information about the nuclear radii. To this end, we consider the HFS for 
the s-electronic states. In this case, first, the effect is maximal, and second, the 
quadrupole effect of the nuclear deformation is absent. 


2 Physical premises 


In the case of the Ml transitions, in the long-wavelength limit, the electromagnetic 
interaction of the nuclei are determined by the single formfactor, or the reduced 
current of the nuclear transition J(R). We define it according to |i5j : 


Jo(R)TS°£(i?)) dn R = 


y/2Ii + 1 


iJ(R). 


(1) 
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with being three basic unit vectors, L — the multipole order of the transition, 
and m — the corresponding magnetic quantum number. Below, we keep the name 
of the transition nuclear current for its radial component J(R). This can be related 
to the radiative amplitude of the nuclear transition. Expressing the latter as follows 


F,, = -i J J^(R)A*(R)d 3 R , 


(3) 


and substituting the expression for the vector potential in the case of transitions of 
the magnetic type PUCE] 

A (iM) = -2yfcj L (uR)Tfy(R), (4) 

one arrives at the equation relating the transition nuclear current with the reduced 
transition magnetic momentum (or amplitude) of the nuclear transition [T7] 

/•OO 

{RJ(R)) = / J(R)R 3 dR = y/{L + l)/L{I 2 \\M(ML)\\h ). (5) 
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Now, using in ([5|) the definition of the nuclear magnetic moment ju m 


eh 

2M p c M = 

one arrives at the expression 


(RJ(R)) 
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(7) 


Relations (JS]) (J7|) show that RJ(R) can be treated as the radial distribution of 

the magnetic currents over the nucleus. Therefore, baring in mind application for 
studying nuclear structure, one can advance the past expression to 

(BP+J(R)) = \J^{2I + l)(I + l)/I^-(^) j , (8) 


where M p is the proton mass, {fi)j = (Rj) J — the j-th moment of the radial 
distribution of magnetization over the nuclear volume. 

Energy shift of the state with the total angular momentum F and its projection 
M, nuclear and atomic momenta / and j, respectively, can be viewed as an elastic 
amplitude of the Ml internal conversion (IC) transition 1121 [IT , 10]. The amplitude 
is defined in terms of the interaction of the nuclear and electronic transition currents: 


F c = 


(J(R)s(r)) V(\r - R|) d 3 r d 3 R , 


(9) 


with the photon propagator 

OO 

£>(|r - R|) = e |r ~ R| /|r - R| = , (10) 

1=0 

and the electron current 

s(r) = e'tp'y'ip = 4eKj(r). (11) 

In (1101) . r< (r>) designates the smaller (larger) of the r or R , and r stands for r/r. In 
turn, j L (uR ), hj^uR) are the spherical Bessel and Hankel functions, respectively. 
In the case of a point-like nucleus, r< = r, r > = R. These relations define the 
non-penetration (NP) model in the internal conversion theory. Within the NP 
model, interaction ([9]) with (1101) is fully factorized with respect to the nuclear and 
electronic variables. The whole transition amplitude factorizes into the radiative 
nuclear amplitude and the remaining factor, independent of the nuclear variables. 
The latter just defines the ICC pa ns]. With the account of the finite nuclear 
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size, such a factorization may only be achieved within the framework of the nuclear 
models, among which the most widely used became the models of surface (SC) and 
volume (VC) transition nuclear currents. On the physical ground, the SC model is 
more justified by the Pauli principle, which to a greater extent prohibits motion of 
the internal nucleons. For this reason, the SC model is the cornerstone of known 
tables [181 . It recommended itself as working well in the internal conversion theory. 
Difference of experimental ICC from their tabular values is called as anomalies. The 
anomalies are thus related to the effects of penetration of electrons into the nuclear 
area. We see that the values of thus determined anomalies depend on the nuclear 
model. 

The diagonal matrix element of the IC interaction Hamiltonian H c (Lfi) defines 
the HFS. H c (Lfj,) is spherical tensor of rank L, L, // being the multipole order and 
its magnetic quantum number. In order to maintain succession with the internal 
conversion theory and previous papers mmm, we will consider formulae in 
general case, and in the end we will pass to the mathematical limit. 

Generally, amplitude @ can be considered as a matrix element of H c {L^l)\ 

F c = (b\H c (L/j,)\a ), (12) 

where |a), |f>) are the are the wavefunctions of the initial and final states of the 
atom. In the IC theory, they are usually characterized with the quantum numbers 
| IMjrn) of the nuclear and electronic spins and their projections on the quantization 
axis. On the other hand, the hyperfine shift Spjj is determined by the diagonal 
matrix element in the | FMIj) representation, with F and M being the total angular 
momentum of the atom and its projection. By making use of the Wigner-Eckart 
theorem, after straightforward algebra, the energy shift can be represented as follows 
0 : 

Spij = ^(FM/jltfeOL^I FMIj) = {Iljj] LF}(Ij\\H c (L)\\Ij) . 

The reduced two-bar matrix element in eq. m is defined as follows: 

{I 2 M 2 j 2 m 2 \H c (LM)\I l M 1 j 1 m 1 ) = 

i> 

Recoupling the angular moments in ©, m, one can separate the angular 
variables [15], 16] 

„iw|r—R| . * \ 

— (J(R)j(r)) = 4vHcu £ (jT^*(R)) 

' ' LXM (16) 

x (j T LAf( f )) JL+A(wr<)/i L+A (u;r > ). 


C{I 2 M 2 LM\I 1 M 1 )C(j 1 m 1 LM\j 2 m 2 ) 
\J (2/i + 1)(2 j 2 + 1) 


(I 2 j 2 \\H c (L)\\hj 


(13) 

(14) 

(15) 
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In the case of transitions of the magnetic type, the terms in (1161) with A = 0 are 
relevant. Using definition (jT]) and integrating over the electronic variables, as shown 
in H3 ng, one arrives at the following expression: 

mnc(mij) = f-s- L -iw{i'jir, U}g, <i7) 

p OO 

G= J(R)j(r)j L (ujr < )h L (ur > ) R 2 dR r 2 dr. (18) 

Jo 

Inserting (1181) into (1131) . one arrives at the following expression for the hyperhne 
splitting w of the states with F = I + j and F = I — j: 


w = —4 eujK 


6 nj(2I + 1) 

I(/ + l)(j + l) • 


(19) 


The fact that the hyperhne shift is generically related to the conversion am¬ 
plitude is manifested in the limiting relation between the formulae for the ICC at 
oj —0 and for the HFS [(L2, DJ 


Spij = (-1 y +I ~ F W{IIjj-LF) 


m\— 3 


1 + m) 


(2L + 1)!! 

x lim Ju 2L+L a a (Ml)(2j + 1)(2 L + 1 )(L + 1 )/Lx 

x (I\\M(M1)\\I) 


( 20 ) 


with L = 1, where m and M are the electron and nuclear masses, respectively. In 
(1201) . olq{tI ) is the analogue of traditional ICC a(rl) extended to the case where the 
conversion electron occupies a discrete electronic state D3ED]. It has dimension of 
energy. Note that ea. (l20l) is exact in the sense that its dependence on the nuclear 
model can be fully related to the way of calculation of ICC. In the SC model, the 
magnetic current reads as follows P21 EE]: 

Jsc(R) = D5(R - Ro ), (21) 


Ro being radius of the magnetic current. With Rq = 0, this model includes the NP 
model. The constant D as determined from the normalization condition (ED reads 
as follows: 

+ + (22) 

Then eq. (TTH]) is factorized, resulting in 


w = 


2 ( 2 / + 1 ) 
I{j + 1 ) 


euj 2 K,- 


eh 

~2MpC 
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where k is the relativistic quantum number, and the radial electronic matrix element 


R = j{r)X L {u)R)r 2 dr. (24) 

Jo 

Eq. (|23p can also be obtained from (12011 . In fl24[) . the electronic transition current 

j(r) = g(r)f(r). (25) 

In ()25[) . g(r), f(r) are the large and small radial Dirac wavefunctions of the electron. 
Furthermore, Xl(ujR) is the potential of the electronic transition. Xj J = hi,(ujR) 
within the NP model. In the SC model it reads 


X L {u)R) = [h L (uR 0 )/j L (coR 0 )\j L (ujR). 


The related ICC in eq. (1201) reads as follows: 


Q( L ) = -4k 


a d (rL) = \Q^F \ 2 
a™ 1 1 

L(LTT) C{J - 2 L0|J - 2 


(26) 


(27) 

(28) 


3 The formalism 


In any realistic nuclear model, after separation of the angular variables, the elec¬ 
tronic and nuclear variables are mixed in the radial conversion matrix element (11811 . 
For the purpose of separation of the penetration effects, let us express the nuclear 
transition current in general form as 


J(R) = Jsc(R) + [J(R ) - Jsc(R)} = Jsc(R) + J P {R) (29) 

and, correspondingly, the matrix element (118(1 — in the form 

G = Gsc + G p . (30) 

Here Gsc is the factorized SC model contribution, and G p bares information on the 
penetration effects. It reads as 


G p = 


POO POO 

/ J p (R)R 2 dR[j L (uR) j(r)h L (ur)r 2 dr + 
Jo Jr 

rR roo 

+h L (uR) / j(r)j[Jujr)r 2 dr\ = / J p (R)Q(R) R 2 dR, 
Jo Jo 

rR 

Q{R) = / [h L (uR)j L (ur) - j L (u}R)h L (ur)] j(r)r 2 dr. 
Jo 


(31) 

(32) 
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In the nuclear area, the electronic wavefunctions are well represented by the 
Taylor series mm- 

g(r) = do + a 2 r 2 + ..., f(r) = b\r + b^r 3 + ... for the s states , (33) 

g(r) = air + a%r 3 + ..., f(r) = 60 + & 2?’ 2 + • • • for the pi states . (34) 

2 

By making use of (1331) . (1341) and asymptotic expansions for the spherical Bessel 
and Hankel functions, the electronic current in (1311) can be put down as follows: 

j(r) = ar + c 3 r 3 + ... . (35) 


As a result, we arrive at the series expansion for Q(R ) and, finally, 

,2 


G = 


1 


uj \ 2ir 


(2 1 + 1)(I + l)/^ 0 ,r A 4 


2 M p c 


u 

T 


r+ Y, -{1^3) W - W)‘) • ( 36 > 


i=2,4,... 


Therefore, by comparing experimental HFS with the theory for an atomic level, we 
obtain the root mean square radius R 2 of the nuclear magnetization, and can unam¬ 
biguously predict the HFS value for all the other s- and pi levels. This parameter 

2 

can associate with the effective nuclear radius which is usually extracted from the 
analysis of the data concerning the nuclear hyperfine anomalies. Furthermore, com¬ 
paring experiment with theory for two levels, one can also find the fourth moment of 
the nuclear magnetization distribution and, correspondingly, predict more precisely 
the HFS value for the other levels, etc. Analysis of n levels provides us with the 
even nuclear magnetization moments up to the 2 n-th one, inclusively. 


4 Results 

Atomic calculations were performed by means of the package of computer codes 
RAINE [18] . Fermi nuclear charge distribution was supposed, with typical parame¬ 
ters. Higher order QED corrections were taken into account. A typical value of Ro 
= 7.121 fm was adopted for the radius of the nuclear magnetic current. The vacuum 
polarization potential and the electron selfenergy correction were allowed for as sug¬ 
gested in refs. [23J [23] ■ As a result, for the Is state, the value was obtained Wh = 
5.107 eV, in coincidence with the previous calculation m- Theoretical calculation 
should then be compared to the experimental value. For a long time, the value of 
5.0841g eV was accepted [25] • A somewhat lower value of 5.0863n eV [M], though 
outside the range of the previous error bars, was obtained recently. Furthermore, 
the wavefunctions, calculated in this way, were used in order to fit the experimental 
data. To this end, we incorporated into the calculation the latest QED correction 







[27(;28j A iuqed = —0.0268 eV. Being added to the basic “Coulomb” value of w c = 
5.0825 eV of the bare nucleus, calculated by us taking into account pure electrostatic 
electron interaction with the nucleus, this results in the calculated HFS value of wth 
= 5.0558 eV, which is by 0.6% lower that the experimental value [26]. Then, dealing 
with the difference between the theory and the latter experimental value as shown 
previously, we achieve consensus with the experimental value for this level with the 
root-mean-square radius of the nuclear magnetization R\ s = 6.207 fin. This is the 
only parameter which we can determine by comparison of theory with experiment 
for one level. We dwell on these details in order to show a possible range of scatter 
of the parameters and sensitivity of the model to them. 

At the time being, there are available data concerning the HFS values for the 
upper states. The latest one obtained in the Experimental Storage Ring at GSI is 
0.79750(18) eV [26] . This is noticeably lower in comparison with the earlier mea¬ 
surement at the Lawrence Livermore National Laboratory [29] of 0.820T0.026 eV. 
Regarding theory, our calculation for a bare nucleus, with the Fermi charge distri¬ 
bution with usual parameters, and Rq = 7.121 fin yielded in the Wq = 0.826 eV. 
QED correction may be introduced according to [ 8 j as AEqed = —0.005 eV. The 
most essential difference with the analysis of the Is case arises in the necessity of 
account of the electron-electron interaction A E\ s _ e . In the zeroth approximation, 
we allowed for that through the self-consistent Dirac-Fock method |T 8 |. As a result, 
we arrived at the value of A E\ s _ e = —0.038 eV. This a little exceeds the calculation 
| 8 j, where all the terms up to the second order in the aZ perturbation series were 
taken into account, and the electron screening correction of AE^ e = —0.030 eV 
was obtained If we rely on the latter value, we arrive at the HFS of = 0.791 
eV. Fit to the experimental value of the HFS for this one level only, results then 
in the root-mean-square radius of the nuclear currents of R% s = 6.275 fin. We see 
that R\ s % i? 2 S ; as we neglect the remaining terms in eq. (1361) including R 4 and 
higher moments. We then found that with account of the two radii, R 2 and R 4 , 
a better fitting the multipole moments of the magnetization is achieved if a value 
of Wfa = 0.792 eV is used. Then solving the system of coupled algebraic equations 
(1361) for the Is- (in the H-like ions) and 2s- (in the Li-like ions) states results in the 
reasonably close value of R 2 = 6.118 fm, and R 4 = 6.760 fm. Note that the radii 
obtained satisfy a plausible relation R 2 < R 4 . 

As a test for the self-consistency of the method, we repeated the same calcula¬ 
tions with another value of Rq = 6.2 fm, which is closer to the R 2 obtained and 
therefore, should be more realistic, and provide with a better convergence of the 
series m, as explained above. Indeed, we obtained as the starting values of wj^ = 
5.079 eV, = 0-796 eV, which values are considerably closer to the experimental 
ones. As expected, by solving the system of two coupled equations (|36l) we arrived 
at the same value of R 2 = 6.12 fm, and slightly corrected value of R± = 6.78 fm, 
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Table 1: Representative results of calculation of the i ?2 and R 4 magnetic radii. In the first column, the 
levels used for solving eqs. (l36l) are indicated. With both the levels involved in the analysis, the resulting 
i ?2 and R 4 parameters hold in spite of essentially different starting values of the SC radius Rq. 


level 

Rq, fm 

i? 2 , fni 

R 4 . fm 

Is 

7.12 

6.21 

— 

2 s 

7.12 

6.27 

— 

both 

7.12 

6.12 

6.76 

both 

6.20 

6.12 

6.78 


differing by 0.3%. Now this was achieved already without additional varying the 

w?? value. The scheme of calculation is illustrated in Table 1. 

Let us see how this works for the 2pi state. Direct data concerning the HFS value 

2 

for the upper states, are not available at the time being, to the best of our knowledge. 
Proceeding in the way similar to the 2s case above, we obtained the Coulomb base 
value of w c = 0.28509 eV for the 2 pi state. Account of the R 2 - and AL-terms in eq. 

2 

(13611 lowers this value to 0.28498 eV. Allowing for the correction A w e - e = —27.19 
meV, owing to the electron-electron interaction, and the QED correction A wqed 
= —0.26 meV [30], we arrive at the resulting value of u> 2 P = 0.25753 eV. This may 
be compared to the theoretical value of 257.84(5) meV predicted in Ref. [30] • 


5 The Bohr-Weisskopf effect 

In a series of papers, it was proposed to cancel the BW effect in linear combination 
of the HFS’s of the two levels (specific difference) with a parameter C (0 E] and 
refs, cited therein): 

A'w = w 2s - (w ls , (37) 

with the ( value of 0.16886. However, a model-dependent character of such relations 
was shown above. Eq. (13711 is derived from the conventional Weisskopf model with 
constant nucleon wavefunctions inside the nuclei. Generally, the specific differences 
(13711 are calculated with the same accuracy like each of the terms, that is ~1%. It 
is worthy of noting here that the exact BW value is unobservable experimentally: 
this is a purely estimated value. The ultimate definition of the BW effect can 
be formulated as the difference between the experimental value of HFS and that 
calculated within the NP model: 

b K = (w N p - w exp )/w NP . (38) 
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The latter model would unambiguously and adequately predict the anticipated effect 

for the corresponding point-like nucleus. This model excludes influence of nuclear 

dynamics. At the same time, in this definition static nuclear properties like charge 

distribution influence the calculated NP value, together with QED corrections. We 

will not further separate out these effects in the definition (1381) . We summarize that 

the calculated NP values are 5.162, 0.811, and 0.25828 eV for the HFS values of the 

Is-, 2s- and 2pi levels, respectively. These result in the BW effect of 1.73, 1.47 and 
2 

0.29 percent, respectively. It is these values of the BW effect which predetermines 
the difference between various models. We recall that the static nuclear properties 
contribute at the level of up to 10% to the ICC values in the case of heavy nuclei 
(e.g. [3l| and refs, therein), that is even more than the BW effect (1381) . 

Our above fit to the data on the HFS of the both Is- and 2s levels results in 
the C value of 0.1746, that is 3% higher than that cited above. This discrepancy 
is just what was fairly expected in the light of what was said above. We conclude 
that parameter £ in the specific difference (1371) is neither observable, nor model- 
independent enough, to be predicted up to five digits by making use of only one or 
two experimental HFS values for the element. 


6 Discussion 

Consrcutive theory of HFS, including the description of the anomalies caused by 
the Bohr—Weisskopf effect, is presented previously. That is founded on the clas¬ 
sical theory of anomalous internal conversion. First, the theory shows the kind 
of information which can be obtained from studies of HFS’. This is even multi¬ 
pole moments of the nuclear magnetization distribution. They can be extracted by 
means of solving the system of coupled equations (1361) for several electronic config¬ 
urations, for which the experimental HFS' are available. Analysis of the data on 
the HFS in the H-like and Li-like ions of 209 Bi resulted in the second- and, for the 
first time in internal conversion theory, the forth-momentum radii: i ?2 = 6.12 fm, 
and = 6.78 fm. Dependence of the HFS on the and higher moments of the 
magnetization distribution turns out to be a specific feature of HFS. Traditional 
IC theory assumes that there is only one parameter, R 2 , which can be determined 
from comparison of theory with experiment. The same situation ican be noted in 
the contemporary analysis of the anomalies in the HFS, observed within the iso¬ 
topic chains: usually the only nuclear parameter, equivalent to R 2 , is discussed 
(e.g., IQ E|). Therefore, dependence of HFS on the higher magnetization distri¬ 
bution moments, demonstrated above, is an extension of the internal conversion 
theory. Furthermore, determination of higher radii from experiment is impossible 
in the absence of precision atomic calculations with proper account of higher order 
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QED corrections. Hence, the above way of analysis simultaneously offers a critical 
test of QED, otherwise impossible. The stringency of the test is only limited by the 
number of the electronic configurations of the atom of the same element for which 
data on the HFS are available. Our fit confirmed validity of the QED corrections 
calculated previously, within the framework of the latest experimental data. 

As we saw above, the corrections at the level of ~1 meV to the HFS of the 2s 
level were crucial for treating the BW effect, which comprises ~0.1% of the shift. 
This clarifies the sensitivity of the method to the parameters. The i ?4 value was 
obtained within the accuracy of about 0.3%, and the R 2 value was most stable 
against fitting, not changed at all. This estimation shows more definitely the limits 
of the theory, as well as its abilities and prospects concerning future experimental 
data while they will be available. 

These conclusions are derived by means of consecutive treatment of the BW 
effect, instead of fighting it, e.g. through combining the specific differences. The 
principal defect of the latter method is the loss of information concerning the nuclear 
structure, as this information is just conveyed via the BW effect. Another defect 
is that the recipe itself is erroneous, being founded on a model consideration. This 
means that for more than well two decades the research in this field has been 
actually going in the wrong direction, when fighting the BW effect instead of using 
it fruitfully for retrieving information from the data. Narrowing of the theme is also 
redicing the discussion about the nuclear anomalies to the single parameter - the 
effective nuclear radius, losing sight of R 4 and other moments. 

Prediction is also made for the HFS of the L2-electron, in which case data are 
not available up to date. Future experiments should manifest real relation between 
the discussed factors. 

To summarize, the new nuclear-model independent way of treating the HFS is 
proposed, proceeding from the theory of anomalies which were first observed in 
ICC. The way is based on simultaneous analysis of experimental data concerning 
HFS for several electronic configurations of the atom. The method can be used as a 
stringent test of QED. Consequently, further development of the experimental basic 
research is needed, aimed at measuring the HFS on ions with various few-electronic 
configurations with a small number of electrons, where calculation of the electronic 
interaction is more feasible. For this purpose, the storage rings are suitable which 
are available or being in reconstruction e.g. at the Lawrence Livermore National 
Laboratory, GSI Darmstadt, Lanzhou in China. From the theoretical point of 
view, the QED corrections, nuclear recoil effect are of great importance for treating 
future data. A more consecutive way of incorporation of these effects into the 
calculation of the atomic wavefunctions looks quite feasible and should be pursued 
in further invest ligation. 
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